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Abstract 

We present theoretical study of electronic structure (spectral functions and Fermi surfaces) for 
incommensurate pseudogap and charge density wave (CDW) and commensurate CDW phases of 
quasi two dimensional diselenides 2H-TaSe2 and 2H-NbSe2. Incommensurate pseudogap regime is 
described within the scenario based on short-range order CDW fluctuations, considered within the 
static Gaussian random field model. In contrast e.g. to high-T c cuprates layered dichalcogenides 
have several different CDW scattering vectors and electronic spectrum with two bands at the Fermi 
level. To this end we present theoretical background for the description of multiple scattering 
processes within multiple bands electronic spectrum. Thus obtained theoretical spectral functions 
and Fermi surfaces are compared with recent ARPES experimental data, demonstrating rather 
good qualitative agreement. 

PACS numbers: 71.10.Hf, 71.30.+h, 71.45.Lr 
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INTRODUCTION 



Quasi-two-dimensional dichalcogenides TX 2 (T=Nb,Ta,Mo,Hf; X=S,Se) and their differ- 
ent polymorhic modifications long time ago attracted the attention of scientific community 
This was connected with: (i) early suggestions to look for high-T c superconductivity in 
layered compounds; (ii) the discovery of phase transitions with formation of charge density 
waves (CDW) In particular, in 2H-TaSe 2 (2H means - hexagonal structure with two Ta 
layers in the unitcell) the second order transition into incommensurate CDW phase is ob- 
served at temperature =122. 3K. At 90K there is another transition to commensurate CDW 
phase l|, l2j . In 2H-NbSe2 transition to incommensurate CDW phase happens at much lower 
temperature of 33. 5K and no commensurate CDW phase is observed. 

Above the temperature of incommensurate CDW transition in these systems there might 
be the range of temperatures where short-range order CDW fluctuations with finite, but 
pretty large, correlation length £ may exist due to low- dimensional nature of these systems 
(and in analogy with antiferromagnetic fluctuations in cuprates). This indeed is experimen- 
tally observed in angular resolved X-ray photoemission (ARPES) experiments 3j-[5]. 

In this paper we present band structure calculation results for 2H-TaSe 2 and 2H-NbSe 2 
with analysis of possible topologies of the Fermi surfaces upon doping, showing possibility 
of formation of "bone"- like Fermi sheets. Further we describe the details of theoretical 
description of multiband electronic multiple scattering on CDW in a multiple band systems, 
as applied to pseudogap, incommensurate and commensurate CDW phases for both 2H- 
TaSe 2 and 2H-NbSe 2 . As an outcome, we obtain spectral functions and Fermi surface maps, 
which are compared with a number of recent ARPES results 3|, |4j. 

BAND STRUCTURE 

The 2H-TX 2 layered compounds have hexagonal crystal structure with the space group 
of symmetry P63/HIHIC with lattice parameters for Ta system a=3.436 A and c=12.7 A. 



Corresponding Wyckoff positions are for Ta 2b (0,0,0.25) and Se 4f (±,§,0.118) |2j. For- 
mal electronic configuration of Ta is d 1 . To calculate electronic structure of the compound 



the linearized muffin-tin orbitals method (LMTO) js] with default settings was employed. 
Obtained band structure and Fermi surfaces are in good agreement with similar LDA cal- 
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dilations by other authors [7|. We do not present any LDA results on 2H-NbSe2 since its 



crystal structure [2[ and corresponding band structure are very close to those of 2H-TaSe2. 

In accord with the previous works [7J in our LDA calculations the Fermi level in 2H-TaSe2 
is crossed by two Ta-5d bands with 3z 2 — r 2 symmetry (see Fig. which are well separated 
from other bands. 

Fermi surface (FS) of 2H-TaSe2 has three (in some works two Q]) hole-like cylinders near 
the T-point and two hole cylinders around K-point. Our results are presented in Fig. [lb. 
Here we observe three hole-like cylinders around the T-point. 

During the recent years several ARPES studies detected the experimental FS of 2H-TaSe 2 . 
In particular, in Ref. [s| the electronic structure of the valence band was studied in lT-TaS2 
and 2H-TaSe2. For 2H-TaSe2 it was shown that along T-K direction there are four crossings 
with the FS. Similar picture is also seen in LDA results (Fig. [T^,b). In later ARPES works 
@, i 0] it was observed that FS of 2H-TaSe 2 has more complex topology. Namely, along T-K 
direction there appear the "bone" -like FS sheets. Within the LDA picture one can obtain 
such "bones" by the shift of the Fermi level down by about 0.1 eV (Fig. Ctti,c). 

To improve over simple LDA, in Fig. [2] we show the "experimental" bands with disper- 
sions: 
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with hopping integrals £j obtained from the fit to experimental Fermi surfaces [10j. Cor- 
responding values of ti (in eV) for Ta system are: for the band forming barrels around 
r and K points i =-0.027, ^=0.199, t 2 =0.221, t 3 =0.028, t 4 =0.013, for the band forming 
"bones" t =0.407, ti=0.114, t 2 =0.444, t 3 =-0.033, i 4 =0.011. For Nb system: for the band 
forming smaller cylinders t =0.0003, ti=0.0824, t 2 =0.1667, t 3 =0.0438, t 4 =0.0158, while for 
the band forming larger cylinders t =0.1731, ^=0.1014, t 2 =0.2268, t 3 =0.037, t 4 =-0.0048. 
These bands are used in further calculations below. 
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ELECTRONIC SCATTERING ON CDW 



Commensurate CDW phase 

Consider the schematic picture of the first Brillouin zone for two-dimensional hexagonal 
lattice, shown in Fig. [3j In hexagonal structures under study, the commensurate CDW 
vector is Q = |TM that corresponds to tripling the lattice period. Scattering an electron 
by this commensurate CDW vector returns electron back to equivalent point after triple 
scattering: e(k + 3Q) = e(k). Moreover for hexagonal structures there are in fact six 
equivalent scattering vectors: Qi = (|, Q2 = (f , 7^75) 



7T, Q 3 =1-^*4 3^H, and Q, 
n Refs. pl, llOj, are observed 



—Qi {I = 1, 2, 3). Maxima of Lindhardt function, calculated in Refs. [3|, LlQfl , are observed on 

nn 

these vectors Q. In addition Lindhardt function shows pronounced maxima [3|,110J for vectors 
X = \YK (Xi = (§, 0)tt, X 2 = (§, ^)tt, X 3 = (-§, j= 3 )7c and % = -X t (I = 1, 2, 3)), which 
appear as sums of scattering vectors Q (see Table I of momenta summation). 

Thus an electron with momentum k is scattered by any of thirteen different momenta 
(see Table I): - preserving its initial momentum k; Q (Q f and Qj); X (X/ and X;). 
Thus for one band case to find the diagonal Green's function of an electron G(k, k) and 
twelve off-diagonal (G(k ± Q;,k) and (G(k±X;,k)) single-electron Green's functions we 
have to solve the system of thirteen linear equations ( ITT)) (see Appendix). Such an approach 
can be generalized for a multiple band case with simplifying assumption that intra- 
and interband CDW scattering amplitudes are just the same (see Appendix) Solving these 
equations we can finally find the diagonal Green function G y (k, k) = 1,2- band indices) 
and corresponding spectral function: 

A(E, k) = —Im V G u (k, k) (2) 

i 

determining the effective electron dispersion. 



Incommensurate CDW phase 

As was pointed above at temperature T = 90 K 2H-TaSe2 (and 2H-Nb5'e2 at 33. 5K) 
undergoes phase transition into incommensurate CDW phase with scattering vector Q ~ 
0.58 — 0.6TM. Similar to the commensurate case discussed above, this vector corresponds 
to six independent scattering vectors Qz, Q2 I = 1,2,3. Let us consider single scattering 
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TABLE I: Table of scattering vectors summation 
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of an electron with momenta k near the FS by vector Q(Q;, Qi). For general values of k, 
such scattering act moves an electron quite far away from the FS, the only exception is an 
electron in the vicinity of the "hot-spots" where e(k+Q) = e(k)). Among multiple scattering 
processes most probable will be successive scattering processes by vectors Q; and Q/ since 
in this case the scattered electron is back again to initial point with momenta k close to the 
Fermi surface. To this end further we will work in the so called two-wave approximation, 
when scattering act consists of two successive scattering processes by vectors Q; and Q;. 
Assuming that scattering amplitude is the same for intra- and interband transitions, for 
diagonal Green function's we obtain (corresponding diagrammatic representation see in Fig. 
H: 

G^(k,k) = ^(k)^ + < 7 i (k)E^G^'(k,k), (3) 



where S = A 2 £^(k + Qi) + 0*(k - Qi)) and ^(k) = E _ e ^ )+i5 
function for the n-th band. Summing Eq. (j3]) over i, one can get: 



VG^'fk k) 9 ' {k) 

(k.kj-^^ (k) . 



is bare retarded Green's 



(4) 



Then using Eq. ([3]) again we obtain: 



G-(k,k)=^(k)^ + T ^g^ : 



(5) 



which can provide us with the spectral function ([2]) for the case of incommensurate CDW 
scattering. 
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CDW pseudogap fluctuations 



Above the temperature of incommensurate CDW transition there is no long-range charge 
ordering, but due to low- dimensionality of the system there are rather well developed short- 
range order CDW fluctuations, with finite correlation length £ and characteristic wave- 
vector Q which becomes rather quick commensurate with Q = |TM [4] as temperature 
lowers. In analogy to incommensurate electronic CDW scattering we will employ two-wave 
approximation with pair of vectors (Qz,Qz). Diagrammatically such scattering processes 
are show in Fig. where three types of interaction lines correspond to three characteristic 
transfer momenta I = 1, 2, 3. 

Let us assume fluctuations Gaussian. Then averaging over such fluctuations corresponds 
;o all possible interconnections of incoming and outgoing interaction lines of the same type 



12 



14j], producing appropriate effective interactions, assumed to be of the form discussed 



in these works. For the case of high enough temperatures one can neglect dynamics of 
fluctuations and average over static random field of Gaussian pseudogap fluctuations 



12 



14j. 

Let us mention that the number of different diagrams is defined by product of number 
of ways to interconnect vertices of type 1, 2 and 3. Since only outgoing and incoming 



ines of each type can 



3e connected, combinatorics corresponds to incommensurate case 



12J. Following Refs. 



12 



14| we use the basic property of the diagrams of this model: 



any diagram with crossing interaction lines is equal to some noncrossing diagram of the 
same order. Thus only noncrossing diagrams can be considered, while contributions of all 
diagrams can be accounted by combinatorial prefactors. And for each type of interaction 



mes 

mm 



,2,3) we will have its own incommensurate combinatorial prefactors, same as in Refs. 



Recurrent procedure for Green's function: single band case 



Within some straightforward generalization of the approach of Refs. [12|, [13( , for single 
band case one-electron Green's function can be obtained via recurrent procedure, which is 
shown diagrammatically in Fig. EJ Here n\ is the number of interaction lines of type I 
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surrounding the "bare" electron line. Analytically this procedure can be written as: 
where ni, n 2 , n 3 - are even and 

S ni + l,n 2 ,n3 = ^ s { n l + l)[Cni+l,n 2 ,n3(^ + Ql) + Cni+l,ji 2 ,n3 (k-Qi)]- (7) 

The other self-energies £ in ([6]) can be found similarly to ([7]), but n 2 or n 3 should be increased 
by one and vectors Q 2 or Q 3 should be added (subtracted) to (from) k, while 

G r ni+l 1 n 2 ,n 3 ( k =•= Ql) = fi , n 1 +l,n 2 ,n 3 ( k =•= Ql) ~~ S n 1 +2,n 2 ,n 3 , (8) 



and 



Here 



^711+2, n 2 ,n 3 

A s(m + 2)G„ 1+ 2,n 2 ,n 3 (k). (9) 



2"™ (k) " E-e(k)+inv(ky (10) 
and k = l/£ is the inverse correlation length of pseudogap fluctuations, n = n\ + n 2 + n 3 , 

v(k) = \v x (k) + v y (k)\, v Xi y(k) = |jr— are projections of quasiparticle velocities. 

For the case of incommensurate fluctuations, combinatorial prefactors are: 

, 2±1 for odd k 
s(n) = { 2 (11) 
| for even k 



This recurrent procedure is applied in analogy with refs. [12|, [13|]. As a first step, one 
takes large enough n = n\ + n 2 + n 3 , for example even, and assume that all Green functions 
G nuri2n3 with even n 1; n 2 , n 3 , such that n = n\ + n 2 + n 3 are equal to zero. Then from 
(EJ) one can find that all £ ni>n2)n3 for the same indices are equal to zero too. Then using 
recurrent procedure, one can get all new values for G nitn2>n3 with even m, n 2 , n 3 , such that 
rii + n 2 + n 3 = n — 2, and repeat the recurrent procedure until one obtains physical Green 
function: 

G(k) = GVo(k). (12) 



Multiple bands pseudogap model for quasi two dimensional hexagonal structures. 



In hexagonal diselenides TaSe2 and NbSe2, as we have seen, the Fermi level is crossed 
by two bands. Thus, our recurrent procedure should be generalized for multiple bands. We 
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will follow Ref. 1111. devoted to description of possible pseudogap behavior in iron based 
superconductors, and assume that intra and interband pseudogap scattering amplitudes are 
identical. This simplifies further analysis and the recurrent procedure for diagonal elements 
of the general matrix (over band indices) Green's function G 1 ^ can be drawn diagrammati- 
cally as in Fig. [7J For our two band model each of the band indices run over two possible 
values, and there is summation over all possible values of indices p, m, I in the vertices (Fig. 
[7J. Thus, the self-energy in these diagrams has no dependence on band indices at all and we 
can obtain the recurrent procedure for G ni ,n 2 ,n z = , ;^ni,n 2 ,n 3 > which is identical to Eqs. 
dl)])-© as in single band case, with only replacement: 

<7n™(k) = Xy niiB2 , n = E E-e^^nv^). n = ^ + ^ + *s- (13) 

and at the end of the procedure we define the physical matrix Green function as: 

G i ^(\c) - a 1 (k)S- ■ + ^Ao( k ) S #Uo( k ) (u] 
{k) ~ ^°>o( k ^ + ! _ E^ Q)0> o(k) ' (M) 

where £ = S x + So,i,o + ^o,o,i- This Green's function allows us to find the spectral function 
(T5]) in the presence of CDW pseudogap fluctuations. 

RESULTS AND DISCUSSION 

In our calculations we used rather typical estimate of CDW potential A =0.05 eV, and 
for correlations length of pseudogap fluctuations we assumed the value of £=10a (where a 
is the lattice spacing). To mimic for the experimental ARPES resolution we broadened our 
spectral functions with Lorentzian of width 7 = 0.03eV, practically it means that we made 
substitution E — » E + ry during all calculations. 

In Fig. [8] we show spectral function maps along high symmetry directions with k z =0 
for 2H-TaSe2- Upper panel shows spectral function map for incommensurate pseudogap 
phase obtained within our pseudogap model. In general it reminds bare "experimental" 
dispersions plotted on Fig. [2j However, here we see some additional broadening of initial 
spectra. These broadened regions of spectral functions with lower intensity represent regions 



pseudogap formation. Why do we speak about regions? In contrast to cuprates [15[ where 



we have finite and rather small number of "hot-spots" here we have almost infinite number 
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of "hot-spots" and it is the interplay between all of them, which leads to the formation of 
such regions. But still dispersions here does not have any obvious discontinuities. 

Middle panel of Fig. [H] displays the case of incommensurate CDW phase. Now we see 
that regions previously covered with pseudogap have clear discontinuities - gaps and also 
many shadow bands. When we transfer further to commensurate CDW phase (lower panel 
of FigJH]) those gaps become even stronger and we can see much more pronounced shadow 
bands. 

Figures [9HTT1 show spectral functions maps in the vicinity of the Fermi level along cuts 
shown on Fig. [Tb. In all figures upper row represents experimental data of Ref. 4j, while 
lower shows our theoretical results. Generally speaking for all phases: incommensurate 
pseudogap Fig. EH incommensurate CDW Fig. [TO] and commensurate CDW Fig. HU we 
obtain quite good qualitative agreement of theory and experiment for the number of bands 
crossing the Fermi level, their position and relative intensity. 

In Fig. [12] we present comparison of experimental and theoretical Fermi surfaces for 
2H-TaSe2- In the middle part of Fig. [12] we show experimental ARPES data from Ref. (J). 
Data at 180K corresponds to the pseudogap phase, while those at 30K are for commensurate 
CDW phase. 

Upper panel of Fig. [12] shows our theoretical Fermi surface in the incommensurate pseu- 
dogap regime for 2H-TaSe2- In general it more or less reminds LDA Fermi surface from Fig. 
[lb. However there are obvious signatures of partial destruction of the Fermi surface sheets. 
Namely, cylinder around K-point and "bones" along K-M direction are partially smeared 
out. It is seen that this picture agrees well with the experimental ARPES data of Ref. jsl, S]. 

For the commensurate CDW phase (lower panel of Fig. [T2l Fermi surface stays close to 
that obtained in LDA and shown in of Fig. [lb. In contrast to incommensurate pseudogap 
phase Fermi surface sheets here are more sharp both in experiment and in theory. The 
cylinder around K-point is now continuous. In the middle of the "bones" we observe the 
start of formation of small triangles as shown in the center of middle panel. Thus, here 
in commensurate CDW phase of 2H-TaSe2 we again obtain an overall agreement between 
theory and experiment. 

In Fig. [13] we show comparison between experimental (middle panel) and theoretical 
Fermi surfaces (lower and upper panel) for 2H-NbSe2. Experimental data on the Fermi 
surface are available only for commensurate CDW phase 5]. Thus we can compare these 
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with theoretical picture shown on lower panel of Fig. [T3j In general, both Fermi surfaces 
remind those from Fig. [lb. Account of electron scattering on commensurate CDW leads 
to a small regions of Fermi surface destruction, namely, along T-K and K-M directions. If 
there exists incommensurate pseudogap phase for 2H-NbSe2 at high enough temperatures, 
its Fermi surface will not be changed much by pseudogap fluctuations, as seen in the upper 
panel of Fig. [13]) . 

CONCLUSION 

To conclude, here we presented theoretical results on electronic structure of two- 
dimensional diselenides 2H-TaSe2 and 2H-NbSe2 within different CDW phases. 

First of all we formulated a theoretical approach to account for multiple scattering of 
electrons on different types of CDW, also for the multiple bands case. Further we investi- 
gated spectral functions and Fermi surfaces for the pseudogap, incommensurate CDW and 
commensurate CDW phases. Calculated theoretical spectral functions within the pseudo- 
gap phase demonstrate "hot regions" , where spectral function is additionally broadened. In 
incommensurate CDW and commensurate CDW phases in a place of these "hot regions" we 
obtained opening of the number of gaps at the intersections with rather pronounced "shadow 
bands" . Comparing experimental and theoretical Fermi surfaces in the pseudogap phase we 
observe rather clear signs of partial Fermi surface destruction with formation of a number of 
typical "Fermi arcs" , separated by pseudogap regions. In commensurate CDW phase Fermi 
surfaces are rather similar to initial LDA picture, with pretty small features due to CDW. 
The overall agreement between theory and ARPES experiments is rather satisfactory. 
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Appendix: Scattering on commensurate CDW 



One band scattering 

First let us consider single band case with electronic "bare" spectrum e(k). These "bare" 
electrons are scattered on CDW potential, written as: 

3 

V(r) = 2A^cosQ ; r. (15) 
1=1 

"Bare" retarded Green's function is: 

*< k > = E-M + V (16) 

Let us introduce short notations: g(k) = g\ g(k + Q/) = ff, g(k — Qi) = fj; g(k + X/) = 4>u 
g(k — X;) = (pj. Then for diagonal Green's function G = G(k, k) and twelve off-diagonal 
(F, = G(k + Q i; k); F I =G(k-Q l ,k); $ z = G{k+X h k); $ r = G(k-X h k)) one can get the 
following system of thirteen linear equations (see Table I of scattering vectors summation): 

G = g + gAF 

F 1 = f 1 A(F I + $ 2 + F 2 + G + F 3 + $0 
F 2 = f 2 A{<5> 2 + F- 2 + $ 3 + F 3 + G + F x ) 
F 3 = f 3 A{F 2 + $ 3 + F- 3 + $i + F 1 + G) 
F l = fcA(G + F 3 + $i + F 1 + $2 + F 2 ) 
F 2 = f 2 A(F- 3 + G + Fi + $2 + F 2 + $ 5 ) 
F 3 = f- 3 A{^ + F 1 + G + F 2 + <S> 3 + F 3 ) 
= 4>iAF; $ r = 4>fAF, 

where F = Th=i( f i + F i)- 

Solving Eqs. (TT71) . one can the get diagonal Green's function G = G(k, k): 

r = nK . K = l-a/3-a(/3 + l)-6(a + l) 

y ' l-ap-a(p + i)-b{a + l)-gA[a(p + l)+p(a + l)]' 1 J 

where a = A(f 2 +fr + f- 3 ), {3 = A(f 1 + f 3 + f 2 ), a = A^ l ^ l + <p 2 ) + f 2 (fo + <p 3 ) + h^^)}, 

b = A 2 [/i(0i + 2 ) + / 3 (03 + <k) + M<h + <ti)]- 



(17) 
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Multiband scattering 



Following the approach of Ref. we assume that CDW scattering amplitude A is 



identical for intra and interband transitions. One can define short notations: g % = g t (k) = 
E-ei(k)+i6 > fi(T) = 9 t (k±Q l ); 4> l [{1) = g i (k±X l ), where i - band index. Diagonal and off- 
diagonal Green's functions will have additional band indices The rest of notations are the 
same as for single band case. For diagonal Green's function, in analogy with first equation 
of the system (TPTj) . one can obtain: 

(■ j = g%, + <m J2 X>r + 1-7 r i (19) 

m 1=1 

Introducing G j = Ffa = *f (I) = 9 = Ei^l fiij) = £,#0)5 

4*1(1) = Ei&ltf) anc ^ summing Eq. (fl^l) over i we get: 

3 

Gi = gi + gAj2(Fi + Fi). (20) 
i=i 

The rest of other twelve equations for and $^ are completely equivalent to correspond- 
ing equations of one band case ( II 7p . Thus, we immediately obtain: 

G j = g j K, (21) 



where K is defined in Eq. ( 1181) . However, now the quantities g, /^m, 4>iff\ are summed up 
over all band indices. From Eq. (Tl9|) and ( 120]) . using Eq. (l2Tj) we finally obtain: 

G*' = + g'^—^- = fa + 9 -^(K - 1), (22) 

which allows us to calculate the spectral function (j2J) with the account of scattering on 
commensurate CDW. 
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FIG. 1: LDA bands and Fermi surfaces for 2H-TaSe2- Panel (a) - LDA electronic dispersions. 
Fermi level corresponds to zero. Panel (b) - LDA Fermi surface. Panel (c) - Fermi surface for 
shifted down Fermi level shown on panel (a) with a short line to obtain bone-like Fermi sheets. 
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FIG. 3: Schematic picture of the first Brillouin zone for hexagonal lattice with characteristic CDW 
vectors: Q = |rM (Qi, Q2, Q3) - commensurate CDW vectors. X = \TK (Xi, X2, X3)) - 
vectors after two scattering on Q, which also have significant Lindhardt function maxima. 
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FIG. 4: Diagrammatic representation of diagonal Green function within two- wave approximation 
for electron scattering on CDW. 
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FIG. 5: Example of a diagram with multiple scattering on CDW, where dashed, wavy and zig-zag 
lines incoming lines corresponds to scattering on Qi, Q2, Q3, and corresponding outgoing lines - 

Qi, Q2, Q 3 . 
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FIG. 6: Diagrammatic representation of Green function within single band pseudogap model for two 
dimensional hexagonal systems. (1— >2) denotes the two last terms, where substitutions Q\ — >■ Q2 
and n\ + 1 — > n<i + 1 should be done. 
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FIG. 7: Diagrammatic representation of Green function within multiband pseudogap model for 
two dimensional hexagonal systems. 
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FIG. 8: Spectral functions of 2H-TaSe2- Upper panel - incommensurate pseudogap phase, middle 
panel - incommensurate CDW phase, lower panel - commensurate CDW phase. 





FIG. 9: Spectral functions for incommensurate pseudogap phase along cuts shown on Fig. lc 
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FIG. 10: The same as Fig. [9] but for incommensurate CDW phase. 
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FIG. 11: The same as Fig. [9] but for commensurate CDW phase. 
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FIG. 12: Comparison of experimental and theoretical fermi surfaces for 2H-TaSe2. Upper panel - 
theoretical Fermi surface for pseudogap CDW phase; middle panel - joint picture of experimental 
data pseudogap phase (upper part) and commensurate CDW phase (lower part). Lower panel - 
theoretical Fermi surface for commensurate CDW phase. 
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FIG. 13: Comparison of experimental and theoretical Fermi surfaces for 2H-NbSe2- Upper panel - 
theoretical Fermi surface for pseudogap CDW phase; middle panel - joint picture of experimental 
data LDA Fermi surface (upper part) and commensurate CDW phase (lower part). Lower panel - 
theoretical Fermi surface for commensurate CDW phase. 
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